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Abstract 

The dissipation and decoherence (for example, the effects of noise in quantum compu- 
tations), interaction with thermostat or in general with physical vacuum, measurement and 
many other complicated problems of open quantum systems are a consequence of interac- 
tion of quantum system with the environment. These problems are described mathematically 
in terms of complex probabilistic process (CPP). Particularly, treating the environment as 
a Markovian process we derive an Langevin-Schrodinger type stochastic differential equation 
(SDE) for describing the quantum system interacting with environment. For the ID randomly 
quantum harmonic oscillator (QHO) model L-Sh SDE is a solution in the form of orthogonal 
CPP. On the basis of orthogonal CPP the stochastic density matrix (SDM) method is devel- 
oped and in its framework relaxation processes in the uncountable dimension closed system 
of " QHO+environment" is investigated. With the help of SDM method the thermodynamical 
potentials, like nonequilibrium entropy and the energy of ground state are exactly constructed. 
The dispersion for different operators are calculated. In particular, the expression for uncer- 
tain relations depending on parameter of interaction with environment is obtained. The Weyl 
transformation for stochastic operators is specified. Ground state Winger function is developed 
in detail. 

Introduction 



Recently a great number of papers P concerning the "quantum chaos", i.e. with the quantum 
analogues of classical systems possessing the dynamic chaos features, have been published. The 
investigations are conducted along different directions, such as analysis of distribution of energy 
levels; definition and calculation of quantities, which are responsible for the presence of chaos 
(corresponding to the classical Lyapunov exponents and KS-entropy)in the quantum systems; study 
of localization and derealization of wave functions around the classical orbits; etc. Though in most 
of cases mentioned above one is faced with the necessity to describe a quantum system statistically, 
so far there was not paid much attention to a stochastic behavior of the wave function itself. 

Many problems of great importance in the field of the non-relativistic quantum mechanics, such 
as description of Lamb shift, spontaneous transitions in atoms, etc., remain unsolved due to the 
fact that the concept of environment (which as a rule is random) has not been considered within 
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the framework of the standard quantum mechanics. It is obvious that a quantum object immersed 
into the thermostat (or more fundamentally the physical vacuum) is an open system. Various 
approaches |2] to the description of such systems exist, mainly in application to the problem of 
continuous measurements. One of them is based on the consideration of the wave function as a 
random process, for which a stochastic differential equation (SDE) is derived. But the equation is 
obtained by the method which is extremely difficult for application even in case of comparatively 
simple type of interaction between the system and the environment, so that some new ideas are 
needed Moreover sometimes it becomes necessary to consider the wave function as a random 

complex process even in closed systems (for example, when a classical analogue of the quantum 
system has the features of the dynamical chaos) [HI-IE]- 

To describe the cases mentioned above, recently was proposed a radically new mathematical 
scheme (theory of stochastic quantum processes (TSQP). This scheme is used to describe the non- 
relativistic quantum system, in general case strongly interacting with the randomly environment 
[H]. Note, that the wave function of closed system "quantum object + thermostat" at that time 
is described by L-Sch type SDE. On the base of TSQP was developed stochastic density matrix 
method permitting construction of all the thermodynamic potentials of quantum subsystem. 

In the present paper relaxation processes in a closed uncountable dimension system "quantum 
oscillator and thermostat" are investigated within the framework of non-perturbational method. 

In particular is receipted evolution equation for described excitations into thermostat. Exact 
representations are found for both widening and shift (analogous to the Lamb shift) of the ground 
state energy level of the immersed into the thermostat (physical vacuum) QHO as well as the 
entropy of an ground state is calculated. 



1 Formulation of the problem 

We shall consider the closed system "quantum object + thermostat" within the framework of 
Langevin-Schrodinger type SDE 

idt^stc = H^stc, (1.1) 
where ID evolution operator H is assumed to be quadratic over the space variable: 

In expressions <|l.l|) - (ll.2|) the frequency f2(t) is stochastic function of time. Let them have the 
form: 

n 2 (t) = Ql + V2ef(t), (1.3) 

where Oq = const and f(t) is independent Gaussian stochastic process with the zero mean and 5— 
shaped correlation function: 

<f(t)f(t')>=6(t-t'), </(t)>=0. (1.4) 

Constant e is characterized the power of stochastic force f(t). The equation Hl.ljl have asymptotic 
solution ty as (n\x, t) in the limit of t — > — oo: 

* as (ri|M) = e^ (n+1/2)n °V(«|z), 
*(„!*) = (^]^) 1/2 ^ 0xV2 H n (^), (1.5) 
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where <p(n\x) is the wave function of a stationary oscillator and H n {x) is the Hcrmitian polynomial. 
The formal solution of the problem l|l.lfl - lll.5|) ma y be written down explicitly for arbitrary f2(t). 
It has the following form: 



V stc (x,t\{0) = -L eX p{lfr- 1 a; 2 }x(^r), (1.6) 

where the function x(Di t ) satisfies the Schrodinger equation for a harmonic oscillator on the 
stochastic space-time {y, r} continuum: 

l ^ = -2W + ~r x ' (L7) 

where 

t 

V=-, e(t)=r(t)e*rW, r = ^, 7 (t) = ^ / *L 

— oo 

The function £(t) is defined from the classical homogenous equation of motion for the oscillator 
with the frequency il(t) 

£+n 2 (t)Z = Q. (1.8) 

Taking into account (|1.6|) and well known solution of autonomous quantum harmonic oscillator 
[TU| for stochastic complex processes which are described closed strong interacting system " quantum 
object+environment" we can write following expression: 



* stc {m\x,t; {£}) = 



i r^ i/2 



2™n! y Ttr(t) 



t 

exp 



The solution of (|1.9|) is a random complex process defined on the extended space 3 = R 1 ® 
where i? 1 is one dimensional euclidian space and R{£} is the corresponding functional space. It is 
easy to show that the mentioned complex random processes are orthogonal. Taking integral over 
the space R 1 we get: 

*stMx, t; {t})V* stc (m\x, t; {£}) dx = 5 nm , (1.10) 



where the symbol * means the complex conjugation. The relation (|1.1U|) shows that the closed 
odd dimensional system " quantum oscillator+environmcnt" is described in terms of the full orthog- 
onal basis of quadratically integrable functionals of the space L 2 . The last fact is very important 
for the further strong mathematical constructions for the statistical parameters of the system. 
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2 Stochastic density matrix method 



The quantum system is impossible to isolate from the environment. This is a principal problem 
if taking into account that even for the ideal isolation nevertheless any system is located in the 
fundamental physical vacuum. Remember that many important features of atomic systems, like 
Lamb shift of energy levels, spontaneous transitions, etc., are explained by vacuum fluctuations. 
So, the processes in quantum systems are to some extent irreversible. For the investigation of 
irreversible processes the non-stationary density matrix representation based on quantum Liouvile 
equation [TI| is often used. However the application of this representation has restrictions. It 
is used for the cases when the system before the interaction switch on was on thermodynamic 
equilibrium state and after the interaction was applied, its evolution is adiabatic. Below in the 
frames of considered model the new approach is used for the investigation of statistical properties 
of irreversible quantum system without any restriction on the quantities and the rate of interaction 
change. This method is based on bilinear form is constructed from orthogonal complex random 
processes. Below it is referred as the stochastic density matrix method. 

Definition 2.1. The stochastic density matrix is defined by the expression: 

oo 

p stc (x, t; t'; {?}) = 4 m) /&, t; {£}!*', f ; {£'}), (2.1) 

m=0 

where p^J(x, t; {£,}\x' , t'\ {£'}) is partial stochastic density matrix and defined by help of bilinear 
form: 

P { £>(xMZ}W >W}) = *Mm\x,t;{Z})** stc (mWA{e})- (2.2) 

In the expression of Wq has the meaning of the initial distribution over quantum states 

with energies E m = (m + l/2)r2o, until the moment when the generator of random excitations 
is activated. Integrating }2.1\) over euclidian space and taking into account M.1U\) . we obtain the 
normalization condition for the weight functions: 

oo 

X>£ m) =l, ^ m) >0. (2.3) 

m=0 

Below we'll define the mean values of various operators. Note that when averaging on extended 
space S the order of the integration is important. If the integral is taken first on R 1 space then 
on R{£}, the stochastic density matrix becomes equal to unity. This means that in the extended 
space all conservations laws are valid, in other words the stochastic matrix in this space is unitary. 
Else if we take the integration in the inverse order, we get another picture. After the integration 
on R{£} the obtained density matrix describes quantum processes in Euclidean space R 1 . Its trace 
is not unity, in general. This means that the conservation laws are not valid already. This can 
be explained by the fact that the system has been in a strongly non-equilibrium state and after 
relaxation its parameters have been significantly changed. 

Below we'll be interested in quantum subsystem processes, hence the integration first on R{£} 
and then on R 1 is supposed. 

Definition 2.2. The expected value of the operator A(x,t\{£}) in quantum state with the index m 



4 



is: 



A,. 



lim 



1 



Sp x Sp^Ap^J \, N m (t) = Sp x Sp {e} p% 



(to) 



(2.4) 



t^ + oo{N rn (t) 

The mean value of the operator A(x,t\{t;}) over the whole ensemble of states will respectively be 
given by: 



A = lim 

t — >+oo 



N(t) = Sp x Sp^yp s tc 



The operation Sp^ in (|2.3|l and (|2.4|l is denned by functional integral representation 

Sp m {K( Xl t;{^}\x',t';{a} = \[^- J K(x,t; {£}\x' ,t; {£}) D{£}, 
and correspondingly the operation Sp x is defined as a simple integration 



Sp x {K(x,t;{Z}\x',t';{?}} = J^ / K(x, t; {d}\x, t; {£'}) dx 



(2.5) 



(2.6) 



(2.7) 



If one wishes to have the quantity describing irreversible behavior of the system, it is necessary to 
change definition of entropy. 

Definition 2.3. The von Neumann entropy, the standard measure of randomness of statistical 
ensemble described by density matrix , which is defined as: 



S N (e,t) = -j^-Sp x {plnpj 



(2.8) 



where N(t) = Sp x p an d p = Sp^j{p s t c }- 

The definition \2. of entropy is correct for quantum information theory and also is agree with 
the Shannon entropy in the classical limit. 

Often it is interesting to know entropies of isolated quantum state (partial entropy) 



S 



(to) 
N 



fe*>=-]^M> (m, *> (m) }' 



(2.9) 



where N m (t) = S Px p^ and p( m ) = 



Definition 2.4. The entropies can be defined another forms. The total entropy may be calculated 
by formula 



S G (e,t) 



N(t) 



Sp x \ Sp [i} 



and the partial entropy correspondingly by formula 



e,t) 



1 



N m (t) 



Pstc In p 8 tc 



Sp x < Sp {i} p\ t J\np 



(2.10) 



(2.11) 



Before proceeding further to the calculations of the physical parameters let us write down the 
general form of the partial stochastic density matrix: 



o 



irr{t)r{t') 



— exp< % 



i(m+\) 



fit 



dp. 



dp, 



■5 



r(t') * . 



'(f) 



fio - 777 )i?n 



'fin 



-(*'). 



(2.12) 



3 Calculation of thermodynamic potentials 



Now let us turn to the calculation of the ground state equilibrium entropy, which is defined by the 
expression (|2.11|) : 

?(°)m i;,„ <*"••( \ / 1 (3.1) 



5^ 0) (A) =Jim 5W(A,t). 



i — >+oo 



The solution of equation 1)1.8(1 may be presented to the form: 

( £o(t) = exp(ifint), i < in = -oo, 

£o(*n)exp{/$(*')dt'}, t>t 0) 



£(*) 



(3.2) 



where $(<) is some complex function. 

After the substitution (|3.2|1 in the (|1.8(l we can define the following nonlinear SDE for a function 

$ + $ 2 + fi 2 + y/2ef{t) = 0, $(to) = &(t )/&(*o) = tfio, (3.3) 

where $ = d t $>. The second equation in the (|3.2|) expresses a condition which guarantees continuity 
of the function £(t) and its first derivative at the t = to. The function <I>(i) is described a complex- 
valued random process due to the initial condition. As a result the SDE i|3.3|) is equivalent to a 
set of two SDE for real- valued random processes. Namely, introducing real and imaginary parts of 

$(t) = Ui(t) + Mi 2 (t), 

we finally obtain the following set of SDE for the components of random vector process u = 



(3.4) 



ui = -u\ + u\ - fig - v^i/(*), f «i(*o) = fle[6(<o)/£o(*o)] = «. 
u 2 = -2uiu 2 , 1 u 2 (*o) = /"i[£o(*o)/fo(*o)] = fio- 

The pair of random processes (wi,m 2 ) are not independent, because their evolution is influenced 
by the common random force f(t). This means that the joint probability distribution: 

P (u,t\u ,t ) = ( JJ £(«;(<) - Uoi)\, u 0i = Ui(t Q ), 



is a non- factorable function. Proceeding from the known evolution equations (|3.4|l . we obtain by 
the standard method the Fokker-Planck equation for Pn (see for ex. ^2] or |13Q which has the 
form: 

<U '" L Q P , (3.5) 



Of 



L (U!,U 2 ) 



+ («? 



d d 

^o) » 1" 2uiU 2 - — 



4ui 



(3.6) 
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with the initial and border conditions: 



P Q (u 1 ,u 2 ;t) 



t=t 



5(ux - u i)S(u 2 - u 02 ), P (ui,u 2 ;t) 



0. 



(3.7) 



Proceeding from 12.8fl . i|2.9[) and i|2.1U[l . we can write down the expression for nonequilibrium 
entropy of ground state: 



S° (e t) - 1 I Na ' a(t) 
^t)- 2 + ^ a(t) 



a=0 



where 



N a (t) = Sp m [l a (t;{Q)}, N a , a (t) = d a Sp U} {l a (t; {£})}. 
In H3.9[l the following notation have been made: 



(3.8) 
(3.9) 



/ a (t;{e}) = ^@exp(-(n + D / „,(/'),//' 
\Ju 2 {t) 



0(u 2 ) 



1, u 2 > 0, 
0, u 2 < 0, 



Now we are ready to calculate the functional traces in expression (|3.9|) . Using the distribution Pq 
it is easy to construct continuum measures in the expressions for iVo(t) = N a (t)\ a —o and N a . a (t), 
which define the entropy of non-equilibrium quantum systems. These integrals may be calculated 
with the use of the generalized Feynman-Kac formula (see for example [§])• So as give the following 
representation: 



+00 +00 



N a (t) 



0(u 2 )Q a (ui,u 2 ;t)duidu 2 



where the function Q a (u\ ,u 2l t) is a solution of the equation: 

dtQ a = L Q a - (a + l)u\Q c 
which is satisfied the following initial and border conditions: 



(3.10) 



(3.11) 



t=t 



5(ui - u i)5(u 2 - U02), Qaiui.U^t) 



I \ u\ I — » + OG 



0. 



Remember that if in the expression H3.10f) the substitution a — is put, the normalization constant 
No(t) will be obtained. If we can calculate the quantity Q a {u\ 1 u 2 ; t) then obviously we can be 
able calculate the function D a {u\ 1 u 2l t) = d a Q a {ui,u 2l t). It is easy to obtain the equation for 
the latter by differentiae the equation (|3.11() with respect to a: 



d t D a = L D a - (a + l)uiD a - uiQ a , 
correspondingly with conditions: 



(3.12) 



D a (ui,u 2 ;t) 



0, D a (u u u 2 ;t) 



0. 
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Introducing the designations Do(ui,u 2 ;t) = D a (u\,u 2 ;t) 



a=0 



we obtain the representation: 



+00 +00 



N 0;0 (t) =N a - a (t) 



9(u 2 )D (ui,u 2 ] t)du\du 2 



(3.13) 



Now the expression for the equilibrium entropy can be obtained: 



4 0) (A) 



1 ^ofo(A) 

2 + <(A) ' 



(3.14) 



where 



jV st (A) = lim JV (t), K (X) = Um JV 0i o(*), A = (n /e^ 3 ) 2 . 



Remember that at the limit t — > +00 the solutions Q Q (ui,U2;i) and Dq(ui, u 2 ; t) turn to their 
stationary limits {111,11%) and Dq(u\,U2) correspondingly. 

Note that the von Neumann entropy (|2.8(l coincides with the entropy expressions in (|2.10() when 
the interaction with the environment is weak, i.e. then f«l. 

To have in view the expressions (|2.5(l and (|2.12(1 by a little transformation for the average 
energy it can be defined: 



E<®(\) = , lim 



1 



-Sp x Sp {6} Hop^tc 



y N (t) 

where operator Hq is determined from l|1.2|) by averaging of H over all: 
Substituting (|3.16(l in (|3.15(l after simple algebra we obtain for the ground state energy: 



with the designations: 



K(X) 



-OO +OO 



iV(f(A) 



>U 2 



-1 



u\ + u\ + X 



2y/X: 



"2 



~}Q t (X;ui,u 2 )dui du 2 



(3.15) 



(3.16) 



(3.17) 



N*(X) 



iti 



"2 



j=Q s (X;ui,u 2 )duidu 2l ui = —-p-, u 2 - , . 
Iu 2 e L ' J e 1 /- 3 



-00 



In the expression 13.18f) the function Qa(ui,u 2 ;t) is a solution of the equation: 
dQo(ui,u 2 ;t) 



dt 



d 2 , - -2 .^9 d 



+ {u x - u 2 + A)-^— + 2u!U 2 -^- + 3ui }Qo{u!,u 2 ; t) 



dui 



(3.18) 



(3.19) 



in the limit of stationary possesses Qq (A; u\, u 2 ) = lim t _^ +oc Qo(ui, u 2 ]t). 
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The energy of n-th quantum state in the thermodynamic limit is calculated similarly: 

4?lW = (« + 1/2) (1 + K(X))n . (3.20) 
As evident from expression (|3.20l) after the relaxation all energetic levels are equidistant. 

4 Uncertain relations, Weyl transformation and Wigner func- 
tion for the ground state 

According to the uncertainty relation, in quantum system the coordinates and momentums can't 
have arbitrary small dispersions. This principle is experimentally verified many times. However 
at the present time as a result of recent quantum technology development the necessity in the 
overcoming of this fundamental restriction and taking the control over the uncertainty relation 
arises. 

The dispersion of the operator Ai is determined by: 



AA(t) = {sp x {pAf)-[Sp x (pA^ 



1/2 



(4.1) 



Now using the expression jUJ it is easy to calculate the dispersions for the operator Ai in the 
extended space S = R 1 ® R^y at the time point t: 



1/2 



(4.2) 



Now, with help of expression l|4.2(l we can calculate the dispersions for the coordinate x and 
momentum p : 

1/2 



Ax(t) 



Ap(t) 



1 



2N (t) 



1 



-3/2 



-oo 
+oo +oo 



2N (t) 



-3/2 



} (ui,u 2 ]t) dui du 2 



'-Qo(ui,u 2 ; t) dui du2 



-oo 



1/2 



(4.3) 



(4.4) 



The product of dispersions for the operators at time point to, when the interaction with the 
environment is switch on describes the standard Hisenberg relation \Ax(t)Ap(t)] = 1/2. It is 

interesting to calculate the uncertainty relation for large time values in the equilibrium limit. By 
averaging t — > oo we obtain: 



Ax st Ap st = lim Ax{t)Ap(t) 
where the following notations have been made: 



1 ^/A X (X)A P (X) 

2 <(A) ' 



(4.5) 



A x (\) = 



-3/2 



<2o'(^; ui,u 2 ) du\ du 2 , 
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A pW = I I ^n^Q^^u^u^d^dui. (4.6) 

-oo * 

where Qg*(A; u-y, U2) is the stationary limit of solution of equation at a = 0. 

It is easy to check that the relations (|4.5[1 - <|4.6[1 differ substantially from the Hisenberg uncer- 
tainty relations. In particular, it allows to control the fundamental relation (|4.5|l using the power 
parameter A, which characterizes the fluctuations of the environment. 

Definition 4.1. We call the expression 

00 

W stc (p,x,t;{0) = 4 m) Wstc(m\p,x,t;{^}), (4.7) 

stochastic Winger function and correspondingly W s t c (iTi\p, x, t; {£}) the partial stochastic Winger 
function. In particular, for the partial stochastic Winger function we get: 



+00 

W stc (m\p,x,t;{£}) = { e^ v ^ stc (m\(x-v/2),t;{^* stc (m\(x + v/2),t;{^})dv, (4. 



Using the stochastic Winger function it is possible to calculate the mean values of the physical 
quantities, which correspond to the operators A: 



+00 +00 

A = a = 



Sp u} {a(p,x,t; {£})pj£(p, x, t; {£})} dpdx, (4.9) 



where stochastic function a(p, x,t; {£}) is defined by Weyl transformation of operator A: 



+00 



a(p, x, *;{£})= / e^y stc (m\(x - v/2),t;{Z})Ay: tc (m\(x + v/2),t;{t;})dv. (4.10) 



Note that the Weyl transformation of quantum operator A after averaging over the random 
process Sp^y{a(p, x, t; {£})} can be used in order to obtain the classical quantity corresponding 
to the quantum operator A. This function is the classical limit H — > (the system of units h = 1 
and c = 1 is used, so the limit h — > actually corresponds to infinite action limit): 



a c i(p,x,t) = lim Sp {i} ^a(p,x,t;{^})Y 
Now we can construct the Winger function for the ground state: 

r0 (x, p; t) = S P{e} { W stc (0\p, x, t; {£})} = 



10 



oo +00 

A f f 1 ( (p — uix) 2 — u\x 2 



2\/-y y ^=exp|-^ ^ ^— |Q(ui,u 2 ;t)dMidu2. (4.11) 

-00 

It is easy to see that function (|4.11|) make sense coordinate x and momentum p simultaneous 
distribution in the phase space at the time t. Particularly with the help of this expression may be 
investigated the relaxation processes and average distribution in the limit of stationary processes. 
Note that such as in the regular case the integration of function W s t c (n\p, x, t; {£}) by phase space 
is normalized per unit: 

+ OC +OC 

W stc {n\p,x,t;{Q)dpdx = l. (4.12) 

Finally it is important to note that for averaging function l|4.11|l it is not to hold expression type 
of gT2>. 



Conclusion 

There are three different reasons which may cause a chaos in the basic quantum mechanical object, 
i.e. the wave function. The first reason refers to measurements performed over a quantum system 
The second reason consists in the more fundamental openness of any quantum system 
resulting from the fact that all the beings are immersed into a physical vacuum 14 . In the 
third place, as it follows from the recent papers [HI EL a chaos may also appear in the wave 
function even in a closed dynamical system. As it is shown in 5], there is a close connection 
between a classical nonintegrability and a chaos in the corresponding quantum system. Many of 
the fundamentally important questions of the quantum physics such as the Lamb shift of energy 
levels, spontaneous transitions between the atom levels, quantum Zeno effect |15j . processes of 
chaos and self-organization in quantum systems, especially those where the phenomena of phase 
transitions type may occur, can be described qualitatively and quantitatively in a rigorous way 
only within the nonperturbative approaches. The Lindblad representation [161 117| for the density 
matrix of the system "quantum object + thermostat" describes a priori the most general situation 
which may appear in the nonrelativistic quantum mechanics. Nevertheless, we need to consider a 
reduced density matrix on a semi-group 0], when investigating a quantum subsystem. This is quite 
an ambiguous procedure and moreover its technical realization is possible only in the framework 
of a particular perturbative scheme. 

A crucially new approach to constructing the quantum mechanics of the closed nonrelativistic 
system "quantum object + thermostat" has been developed recently by the authors of [f3E| from 
the principle of " local correctness of Schrodinger representation" . To put it differently, it has 
been assumed that the evolution of the quantum system is such that it may be described by the 
Schrodinger equation on any small time interval, while the motion as a whole is described by a SDE 
for the wave function. In this case, however, there emerges not a simple problem to find a measure 
for calculating the average values of the physical system parameters. Nevertheless, there exists a 
certain class of models for which all the derivations can be made not applying the perturbation 
theory 
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In the present paper we explore further the possibility of building the nonrelativistic quantum 
mechanics of closed system " quantum object + thermostat" within the framework of the model of 
one-dimensional randomly wandering QHO (with a random frequency and subjected to a random 
external force). Mathematically the problem is formulated in terms of SDE for a complex- valued 
probability process defined on the extended space R 1 ® R{£}- The initial SDE is reduced to the 
Schrodinger equation for an autonomous oscillator defined on a random space-time continuum, 
with the use of a nonlinear transformation and one-dimensional etalon nonlinear equation of the 
Langevin type defined on the functional space R{£}- It is possible to find for any fixed {£} an 
orthonormal basis of complex- valued random functional in the space L,2{R 1 ) of square-integrable 
functions. With the assumption that the random force generator is described by a white noise 
correlator, the Fokker-Planck equation for a conditional probability is found. From the solutions 
of this equation on an infinitely small time interval a measure of the functional space R{£} can be 
constructed. Then by averaging an instantaneous value of the transition probability over the space 
R{£}, the mean value of the transition probability is represented by a functional integral. Using 
the generalized Feynman-Kac theorem, it is possible to reduce the functional integral in the most 
general case, where both frequency and force are random, to a multiple integral of the fundamental 
solution of some parabolic partial differential equation. The qualitative analysis of the parabolic 
equation shows that it may have discontinuous solutions. This is equivalent to the existence of 
phenomena like the phase transitions in the microscopic transition probabilities. In the context of 
the developed approach the representation of the stochastic density matrix is introduced, which 
allows to build a closed scheme for both nonequilibrium and equilibrium thermodynamics. The 
analytic formulas for the ground energy level broadening and shift are obtained, as well as for 
the entropy of the ground quantum state. The important results of the work are calculation of 
expressions for uncertain relations and Wigner function for the quantum subsystem which is strong 
interacting with the environment. 

The further development of the considered formalism in application to exactly solvable many- 
dimensional models may essentially extend our understanding of the quantum world and lead us 
to the new nontrivial discoveries. 
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